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in 

The no-boundary wave function of quantum gravity tends to assign only very small prob 
£ — ■ ability to long periods of inflation. This leads to tension with observations. We study the 

' no-boundary proposal in the context of multi-field inflation to see whether the number of fields 

changes the situation. For a simple model, we find that indeed the no-boundary wave function 
can give higher probability for sufficient inflation, but the number of fields involved has to be 
very high. 



1 Introduction 

The nature of the initial singularity of our universe is an important but unsolved problem in modern 
cosmology. One approach to this problem is to canonically quantize the universe and study its wave 
functions, solutions of the Wheeler-DeWitt equation pQ. There is no universally accepted way to 
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choose, among the many solutions, the wave function that is supposed to describe our universe. 
One way to address this problem is to draw an analogy with the ground state in ordinary quantum 
systems. Hartle and Hawking [2] observed that the latter states can often be written in terms of the 
Euclidean path integral. Thus they suggested using an analogous definition in cosmology, resulting 
in what is called the no-boundary wave function. 

The no-boundary wave function can be interpreted as a superposition of space-time histories. 
It can be approximated by sum-over on-shell solutions, so-called instantons. The sum-over instan- 
tons can have problems of divergence [3] , which may be alleviated by considering complex valued 
instantons [4] , so-called fuzzy instantons [3 [6] . 

As is to be expected in a quantum theory, the no-boundary wave function exhibits classical 
behavior only for a part of the phase space, i.e. only for certain space-time histories. Since we 
find ourselves in a universe which behaves classical on large scales, it is interesting to consider 
those histories in detail. It turns out that for this preferred subset of histories, the no-boundary 
wave-function assigns, to some approximation, a classical probability distribution, called the no- 
boundary measure. Recently, Hartle, Hawking and Hertog investigated the no-boundary measure 
for slow- rolling fields in Einstein gravity 5 [6 ; the present authors studied the no-boundary measure 
for more general examples: scalar-tensor gravity [7] and fast-rolling fields [5]. 

Due to previous studies, it was known that the no-boundary measure disfavors the histories 
which have a large number of e-foldings. If we believe that the early universe had experienced more 
than 60 e-foldings from primordial inflation, the no-boundary measure is not compatible with this 
expectation. 

On these grounds, one may reject the no-boundary wave function in favor of other prescriptions, 
for example Vilenkin's tunneling proposal [9], which may however be regarded as less natural on 
theoretical grounds. One may also try to keep the no-boundary wave function, but 

• try to avoid the need for a long period of inflation such as ekpyrotic universe [TO], string gas 
cosmology [TJJ, or big bounce model [T2"] . 

• justify inflation on anthropic grounds, for example in the context of the multiverse |13) . 

• lift up the inflationary probabilities by considering the volume weighting |14| [15] . 

These proposal each come with their own questions and difficulties attached. Therefore in this 
article we explore another mechanism to lift up the probability for the larger e-foldings in the 
context of the no-boundary measure. We study the no-boundary measure of multiple fields and 
argue that the degeneracy of the large number of fields can lift up the probability of inflationary 
histories without considering any additional factors. Although the Euclidean action does not prefer 
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inflationary histories, there exist many numbers of histories which produce larger e-foldings, if we 
include a number of fields. Then, for large numbers of fields, larger e-foldings can be sufficiently 
preferred. 

In Section (5J we briefly review the no-boundary wave function and explain the reason why 
it does not prefer inflationary histories. In Section [21 the no-boundary measure is extended to 
the multi-field case and N-flation models, and it is demonstrated that the degeneracy can lift up 
the probability of larger e-foldings. We also relate the model with the observations. Finally, we 
summarize the results and discuss possible extensions in Section 3) 



2 No-boundary wave function 

The ground state of the universe, or the no-boundary condition of the wave function of the universe, 
can be described by the Euclidean path integral [2]: 

*IV>X] = J2 [ Vg^V* e- s ^^\ (1) 
r J 

where T denote compact manifold topologies with a single spherical boundary, and the metric g^ v 
and the matter field <f> take the value and x a t the boundary. Here, Se is the Euclidean action: 
we assume Einstein gravity and a minimally coupled scalar field, 



Se = - cfx^+g- —R - -(V0) 2 - V{4>) (2) 



in the units H = G = c = 1. In this paper, for simplicity, we restrict our interests to a minisuper- 
space; for a homogeneous and isotropic universe, the metric becomes 

ds 2 = N(X) 2 dX 2 + a{X) 2 dn 2 3 . (3) 

In this context, it is convenient to view the time integral in the action as a contour integral in the 
complex time plane (for details see for example [5]) and shift the contour in suitable ways. Thus, A 
(which we take to be in [0, 1]) is a parametrization of the chosen contour. The no-boundary wave 
function takes the form 

*[6,X] = J VNVaVcj) e- 3 ^ 11 '^, (4) 

where a(X = 1) = b and </>(A = 1) = \ are the boundary values at the only boundary of space-time, 
and at A — one imposes conditions suitable for a regular metric ("no-boundary conditions"). 

The path integral can be approximated by the method of steepest descent. To lowest order, in 
this approximation we only keep the contribution from saddle points of the exponent: 

$[6, X ]^e^M (5) 
p 



where S^[b, x] is the action for the solution labeled by p with the boundary conditions b, x, of the 
equations 

*="3^ + ™, a = -|a(^ + F W ), (6) 

where ' denotes derivative with respect to r defined by 

r(A) = J X dX'N(X'). (7) 

The on-shell action can be simplified to 

S E = 4tt 2 J dr (a 3 V - J^a) . (8) 

We are especially interested in regions in which the wave function behaves almost classically. In 
the saddle point approximation ([5]), this translates to the condition 

|V5^[6,x]|«|v4 m [6, X ]|, (9) 

which singles out Euclidean saddle points that contribute in the classical region. S^ c and S^ 11 are 
the real and imaginary parts of the Euclidean action and V is a derivation that is defined on the 
superspace. It can be shown [5] that in this region the wave function effectively assigns classical 
probabilities for parts of the phase space, i.e. for some space-time histories. In fact, for a choice of 
contour such that the last part runs along the Lorentzian time axis, the histories that are assigned 
probability are precisely the histories satisfying ([9]), when evolved further in Lorentzian time. Hence 
we will call such histories classical histories. The probability measure for classical histories becomes 

dP[b,x] oc \^[b, X }\ 2 dbd X c± e- 2S v° [b - x] dbd X - (10) 

The no-boundary measure is defined on the boundary of the histories, b and x- However, it is 
convenient to label a history not by the boundary values, but by the initial conditions, for slow-roll 
cases, since the Euclidean action is mostly determined by the initial potential energy: 

where tp = 4>(t — 0). In these cases, one thus has 

dP[p] aex pQ^y) *P- (12) 

as an approximation of the probability to observe a history starting from tp in the Euclidean regime. 
We take the simplest inflation model for single field as one having a quadratic potential [17] 

V(4>) = \r?fi- (13) 
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It is known that except for tp = 0, only tp larger then a critical value has the solution space satisfying 
Equation © 6 . That is, the classical history exists only for 

ip > tp c ~ 0.62. (14) 

Although the no-boundary measure diverges at tp = 0, we will ignore this point. Due to the 
divergence it is impossible to evaluate the derivatives in ©, and it is hard to see how this history 
could be classical. 

In the slow-roll approximation |18) , the total e-folding number obtained from tp is approximated 

by 

N e ~ -8tt / d^> (15) 

where we have ignored a small correction from the location of the end of inflation. This can be 
further reduced for the quadratic potential Equation (|13|) : 

N e ~ 2wtp 2 . (16) 



Then, the no-boundary measure, Equation (|12[). can be expressed as a function of e-folding numbe: 

rfm] ocexp(|-^-)^, (17) 

Here, the Jacobian of the change of variables is ignored because the exponential part is dominant. 
The classical histories exist for 

N e > iV e c ee 2ir{p c ) 2 ~ 2.4. (18) 

Therefore, the no-boundary measure exponentially prefers a small e-folding number (~ 2.4). 

As discussed in Section [TJ there are several possible interpretations of this disagreement between 
the prediction of the no-boundary measure and the inflationary universe. One of the idea was that 
the volume factor should be weighted to measure the probability to observe our past light cone. It 
enhances the likelihood of larger e- foldings by multiplying the volume [5] , 

dP[N e ] a exp (^-^- + 3N e ) dN e . (19) 



V 2 m 2 N ( 

This volume weighting sufficiently lifts up the probability of larger e-foldings, if 

N e > mT 2 . (20) 

In the following, we will study under which conditions the addition of further fields may have a 
similar effect. 



1 Since the expansion during the Euclidean time is negligible, Equation (116ft is still a good approximation. 
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3 No-boundary measure in multi-field inflation 

In this section, we study the no-boundary measure with multiple fields. As simple toy model, we 
take all fields, havingsame quadratic potential type for simplicity as in the earlier works on N- 
flation models 19, 20p where it was shown that any of these fields does not need to take on values 
in excess of the Planck scaljf] and a large number of fields, predicted by string vacuum solutions, is 
needed to get sufficient e-foldings. We take N-flation models as a possible candidate of the multi- 
field inflation models throughout this paper. The basic effect that we will be investigating is that, 
although large amounts of inflation are still suppressed in terms of the action, there are now many 
histories that contribute to a given amount of inflation. This volume in field space competes with 
the suppression due to the action. 

3.1 No-boundary wave function of multiple fields 

The no-boundary wave function in Equation ([1} can be extended to the multi-field case. In the 
minisuperspace of scale factor and a set of scalar fields, it becomes 

V[b,x\ = J VNVaV(f e~ SE[a ^ ] (21) 

where the arrow denotes multiple fields, <p = {4>j}, and a(X = 1) = b and <f>(X = 1) = x an d the 
paramctrization A was chosen such that A = 1 corresponds to the single boundary. 
As discussed in Section [2j we use the method of the steepest descent 

*[6,^~^e- s -M, (22) 
p 

where 

& = -3fjA + ^, a = -fa(£ + V(fl), (23) 

S E = 4^ 2 J dr (a 3 V$ - J^a) , (24) 

where the subscript i denotes the i-th field. 

2 Earlier investigations of the N-flation paradigm 1191 1201 1211 1221 have worked under the assumption that all 

relevant fields are close to their minima and can be described by quadratic potentials. For axions the full potential 

is trigonometric and recently Refs. 1231 showed that the quadratic approximation is unreliable. The reason we follow 

the earlier works rather than the more recent ones is mainly for simplicity. We intend to do further work on for the 

case of more complicated potentials. 

3 Note that, in literature, people use reduced Planck mass scale. 
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3.2 Degeneracy of solution space for the case of a quadratic potential 

We consider a set of uncoupled fields in quadratic potentials with the same mass m% = m [191 120) , 

V{$) = \m 2 $ 2 . (25) 

The equations of motion are 

6 = -3-6 + m 2 $, a = --£-a ( $ 2 + \m 2 $ 2 ) . (26) 

a 6 \ 2 / 

The no-boundary conditions at r = require 6 = 0. Thus we can make the ansatz 

$(t) = 6(r)4>, (0) 2 = 1, (27) 

which requires 6 = (6/6)6. Then the equations of motion for 6(t) and a(r) are precisely given by 
the equations © for a single scalar in a quadratic potential. The no-boundary wave function is 
thus independent of the direction of 6 in field space, 

¥(6, tf) =¥(&,*), X =|4 (28) 

As a result, the classicality conditions in this case reduce to those of the single field case, and one 
can use Equations (|12|) and (|T4]) . Similar to Equation (fT2)) . the no-boundary measure can be defined 
on the initial field space in the slow-roll limit, 

dP[p] cc exp (j-^-s) d<p (29) 

and we have again ignored a Jacobian due to the change of variables, tp = Important physical 
observables, e.g., the e- folding number and power spectrum, do not depend on the angular part of 
if. Therefore, for these observables the effective probability distribution is 

dP[<p] = J d<pdP(0,<p) ~ D(f 2 )exp f|^j^aj ^ (30) 

where |<p| = 1 and D(tp 2 ) is the surface area of the shell \<p\ = ip in field space. 

Now we need to discuss the issue of cutoffs in (l30|) : On the one hand, there are no classical 
solutions below a certain threshold, see (fT5|) . On the other hand, the probability distribution (|3"0"|) 
can not be normalized since it becomes constant for large values of ip, and a cutoff is needed. Here, 
we assume that all fields are below a that is order of the Planck mass ~ 1. This is indeed the 
philosophy of the N-flation scenario [TH1 [201 US H3]- If a field value is larger than the Planck scale, 
it is very difficult to trust their self-consistency of the effective theory; therefore, it is reasonable to 
restrict the field values less than the Planck scale and consider assisted inflation. Furthermore, if 
one assumes the axion fields, the cutoff arises naturally. 
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We will account for the cutoffs by changing D(p> 2 ) appropriately. To estimate D(ip 2 ), we choose 
ifi randomly and count resulting ip 2 . ft can be done by assuming ipi as a uniform random variable 
and using the central limit theorem. We assume that there are Nf number of fields with initial field 
value ifi following 

<PieU(p,a), (31) 

where U(p,q) is the uniform distribution with minimum p and maximum q. For the uniform 
distribution of ipi, the mean and the variance of ip 2 become 



Eyti = -[* 

a Jo 



,2 



dy>i Vt = y, (32) 

a 2 y\\ = E[ri]-E[<fff = l£. (33) 

When Nf is large, a statistical behavior will arise, e.g., it will be more probable to observe ip 2 
around its expectation value. It is described by the central limit theorem which states that the 
distribution of -E[</? 2 ] follows the Gaussian distribution AT (/i, <r 2 ) with the mean /x and the variance 
a 2 , such that 

E\ v 2 ] e N(,,,<, 3 ), (34) 

M . ^, (35) 

= (36) 

This distribution corresponds to the relative number of equivalent choices of <p for a given ip 2 . 
Therefore, 



D(p> 2 ) cc exp I - ^ 2 J> ). (37) 



VA. CJ^p I — 

Then, the no-boundary measure of ip 2 follows from Equation (I30p : 



It represents the Euclidean probability measure and the degeneracy of the configuration space at 
the same time. 

Now we turn our attention to the e-folding number. The degeneracy factor can be described in 
terms of the e-folding number by using Equation (|16[) . 

E[N e ] G N(fi,a 2 ), (39) 
A = 5-^, (40) 




Figure 1: Comparison between jj, and E s i m [N e ] and between a 2 and a 2 im [N e ]. Subscript sim denotes 
the numerical values obtained from 100 simulations for each Nf. We used the uniform random tpi 
following Equation (l3lT) with a = l/y/8n which corresponds to the initial held value below the 
reduced Planck mass. Only Lorentzian evolution is considered in these simulations, since the 
expansion during the Euclidean time is negligible. 



2 16TT 2 a 4 N f . . 

S 45 (41) 

This distribution is confirmed by the simulation in Figure [1] The no-boundary measure also 
can be denoted by the e-folding number: 

) ( 42 ) 
/ 3tt 1 5 . / iV P x J 



= exp - N f 1 dN e . 

y ^2 m 2 N e 8 i \ fi J J 

The degeneracy factor enhances the probability at N e = jx and reduces the probability at N e = N° . 
Therefore, it increases the likelihood of the larger e-folding number. This effect becomes stronger 
as Nf increases and overcomes the probability peak at N£ when 

N f > m' 2 , (43) 

since the parenthesis in the exponential of Equation ([42 p is less than one for N e < ft. Note that 
this condition is comparable to that of the volume factor in Equation (|20p. because when the cutoff 
a is order of the Planck scale, a set of fields satisfying Equation (j43]) will produce enough e-folding 
number by satisfying Equation (|20p . Figure [2] summarizes the main result of this paper. 

3.3 Cosmological viability 

As previously mentioned, we take N-fLation models in the same mass case [THl [20] as possible 
candidates of multi-field inflation and apply the results of the no-boundary measure. 

Using the slow-roll approximation and for quadratic potentials, the number of e-foldings Eq. (|15[) 
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■ exp (m~ 2 - Nf) 




Figure 2: Schematic diagram of the no-boundary measure in Equation (I42[) . For initial fields below 
the Planck scale, the degeneracy factor of the e-folding number has expectation value of the order 
of Nf while its standard deviation is order y/Nf. When Nf > m~ 2 , the peak at N£ which comes 
from the Euclidean probability measure becomes comparable to the peak at jl which comes from 
the degeneracy of the field space. Therefore, the degeneracy factor enhances the likelihood of the 
larger e-folding number. Note that the probability in vertical axis is not normalized and intended 
to represent the relative probability. 

can be rewritten as 

^-2^$. (44) 

We will assume throughout that the observable scales crossed outside the horizon 60 e-foldings 
before the end of inflation N* ~ 60 [26]. 

The power spectrum, the spectral index, the tensor-to-scalar ratio, and the non-gaussianity 
parameteiQ are given by [2U] 

Nf 

t iron ! I 



2 



ns-l - (46) 
r *— ~ — (47) 

5 /NL " i^W^^' (48) 



Here the non-gaussianity is local type. 
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For N* ~ 60 [33], these becomes are given by 

ns-1 ~ 0.033, (49) 

r ~ 0.133, (50) 

^/nl ~ 0(0.01). (51) 
5 

Even though this case does not give significantly large non-gaussianity, but still all these parameters 
are still within viable region of WMAP7 25 . 

We can reduce the formula for the power spectrum in terms of effect single field as 

P K ~ ^m 2 (<n 4 . (52) 

The observed normalization of the spectrum can be taken as P^ 2 ~ 5 x 10~ 5 [25], leading to a 
normalization of 

m ~ 1.3 X lO -6 , (53) 

which is a reasonable mass scale for the quadratic model regardless of the choice of Nr. Therefore, 
from Equation (|43p. the multi-field factor works for Nf > 10 12 which is quite large number of fields 
needed. 



4 Discussions and conclusion 

In this article, we investigated the no-boundary measure with Nf fields. To obtain a very simple 
model, we assumed a quadratic potential with a same mass to. Following the spirit of N-flation, we 
further assumed that the field values of each fields should be less than the Planck scale: a < 1. Then, 
the no-boundary measure effectively has two competing factors (Equation (|38[) ): the Euclidean 
action and the Gaussian degeneracy factor. 

If the number of fields Nf becomes sufficiently large so that Nf ~ to" 2 , then the two factors 
are of similar orders. The peak of the Gaussian factor ip ~ Nt then gives realistic cosmologies 
non/negligible probability. However, the nuber of fields needed to achieve this effect seems to be 
very large. In the simple model we considered, Nf ~ 10 12 number of fields are needed to fit the 
power spectrum of our universe. Since this is a very high number, this seems to suggest that this is 
not an effective way to resolve the tension between no-boundary probabilities and the observations. 
But there are also less pessimistic ways to view the situation: 

1. A large number of scalar fields is not ruled out by experiments if their energy scale is near 
Planck scale. Indeed it has been argued that such a scenario is useful to resolve some problems 
of our phenomcnological universe [16j . Therefore, it may still be a viable idea. 
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2. In our calculation, we assumed a quadratic potential described by a single constant mass, but 
one can certainly question these assumptions. 

(a) The masses can depend on the energy scale or the time, via the running of couplings 
(masses) or the temperature dependence. If the effective mass is sufficiently large when 
the universe begins inflation, and if it decreases as time goes on, then one may explain 
the sufficient e-foldings with reasonable numbers of fields. 

(b) The mass can be different for different fields, e.g. Refs.[20, 22 . If the Euclidean dynamics 
highly depends on the larger mass, then the required number of fields to lift up the larger 
e-foldings can be reduced. This can explain the power spectrum, if the smallest mass is 
order of 10 -6 . 

(c) The quadratic potential should only be regarded as a first approximation of more compli- 
cated potentials(e.g., axions in Refs. [13]). More complicated potentials may change the 
Euclidean cutoffs and reduce the peak of the Euclidean probability, leading to a required 
Nf more in keeping with present models of high energy physics. 

(d) If the fields couple in complicated ways, then again, the Euclidean cutoff structures can 
be changed. Therefore, from the similar arguments, one may reduce the number of fields 
to a reasonable amount. 

To summarize, we have shown that considering the no-boundary wave function for many fields can 
solve the traditional problem of low probabilities for sufficient inflation. But in the simple model 
we have considered, the number of fields is very high, much higher than supported by, for example, 
string theory scenarios. Further research is needed to see whether the high number of necessary 
fields can be brought down in more complicated and realistic models. 

Acknowledgment 

We would like to thank Robert Brandenberger and Ewan Stewart for helpful discussions. DH, DY 
and BHL were supported by the National Research Foundation of Korea(NRF) grant funded by 
the Korea government (MEST) through the Center for Quantum Spacetime(CQUeST) of Sogang 
University with the grant number 2005-0049409. SAK was supported by Basic Science Research 
Program through the National Research Foundation of Korea(NRF) grant funded by the Korea 
government (MEST) with the grant number 2011-0011083. SAK and DY would like to thank the 
organizers of the 2012 Asia Pacific School/ Workshop on Cosmology and Gravitation and the KASI- 
APCTP joint mini-workshop on Science and Technology of BigBoSS. SAK thanks to the hospitality 



12 



of Laila Alabidi, Takahiro Tanaka and YITP, and of CQUcST while this work was being carried 
out. 



References 

[1] B. S. DeWitt, Phys. Rev. 160, 1113 (1967). 

[2] J. B. Hartle and S. W. Hawking, Phys. Rev. D 28, 2960 (1983). 

[3] A. Vilenkin, Phys. Rev. D 50, 2581 (1994) |arXiv:gr-qc/9403010| . 

[4] J. J. Halliwell and J. B. Hartle, Phys. Rev. D 41, 1815 (1990). 

[5] J. B. Hartle, S. W. Hawking and T. Hertog, Phys. Rev. Lett. 100, 201301 (2008) 



|arXiv:071 1.46301 [hep-th]]. 

[6] J. B. Hartle, S. W. Hawking and T. Hertog, Phys. Rev. D 77, 123537 (2008) |arXiv:0803. 16631 
[hep-th]]. 

[7] D. i. Hwang, H. Sahlmann and D. h. Yeom, arXiv:1107.4653 [gr-qc]. 

D. -i. Hwang, B. -H. Lee, H. Sahlmann and D. -h. Yeom. larXlv:1203.0112l [gr-ac] . 
[9] A. Vilenkin, Phys. Rev. D 33, 3560 (1986). 
[10] J. Khoury, B. A. Ovrut, P. J. Steinhardt and N. Turok, Phys. Rev. D 64, 123522 (2001) 



hep-th/0103239|. 

[11] R. H. Brandenberger and C. Vafa, Nucl. Phys. B 316, 391 (1989). 

[12] J. Khoury, B. A. Ovrut, N. Seiberg, P. J. Steinhardt and N. Turok, Phys. Rev. D 65, 086007 
(2002) hep-th/0108187 . 



[13] L. Susskind, hep-th/0302219 



[14] S. W. Hawking and T. Hertog, Phys. Rev. D 66, 123509 (2002) [hep-th/0204212] . 
[15] S. W. Hawking and T. Hertog, Phys. Rev. D 73, 123527 (2006) |hep-th/060209l] . 
[16] G. Dvali, Fortsch. Phys. 58, 528 (2010) jarXiv:0706.2050l [hep-th]]. 
[17] A. D. Linde, Phys. Lett. B129, 177 (1983). 



13 



[18] 



P. J. Steinhardt and M. S. Turner, Phys. Rev. D29, 2162 (1984), D. S. Salopek and J. R. 
Bond, Phys. Rev. D42, 3936 (1990), A. R. Liddle and D. H. Lyth, Phys. Lett. B291, 391 



[19] 



(1992), astro-ph/9208007. 

S. Dimopoulos, S. Kachru, J. McGreevy, and J. Wacker, JCAP 0808, 003 (2008), 



hep-th/0507205 



[20] S. A. Kim and A. R. Liddle, Phys. Rev D74, 023513 (2006), [astro -ph/0605604 
[21] S. A. Kim and A. R. Liddle, Phys. Rev D74, 063522 (2006), |astro-ph/0608186 
[22] R. Easther and L. McAllister, JCAP 0605, 018 (2006), |hep-th/0512102| S. A. Kim and A. 



R. Liddle, Phys. Rev D76, 063515 (2007), arXiv:0707.1982. 



[23] 



S. A. Kim, A. R. Liddle and D. Seery, Phys. Rev L105, 181302 (2010l. larXiv:1005.4410l [astro- 
ph.CO], S. A. Kim, A. R. Liddle and D. Seery, Phys. Rev D85, 023532 f2012). larXrv:1108.2944l 
[astro-ph.CO]. 



[24] A. R. Liddle, A. Mazumdar, and F. E. Schunck, Phys. Rev. D58, 061301(R) (1998), 



astro-ph/9804177 



[25] E. Komatsu et al, Astrophys. J. Supp. 192, 18 (2011). larXiv: 1001 .45381 
[26] A. R. Liddle and S. M. Leach, Phys. Rev. D68, 103503 (2003) |astro-ph/0305263 
[27] L. Alabidi and D. H. Lyth, JCAP 0605, 016 (2006), astro-ph/05 10441 
[28] D. H. Lyth and D. Wands, Phys. Lett. B 524, 5 (2002) [hep-ph/0110002] . 



14 



